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I. Introduction 

The emphasis of this paper is determining necessary and sufficient 
conditions under which the linear partial differential equation 


( 1 ) 


du 

dt 




A .<«> 

, ik <r*x cfy. 

= ! j k 


z b (x) 

.-1 j tfx - 

y- ] J 


= f ( >: * t ) 


can be transformed to become either constant coefficient or of the Kolmogorov 

oo <V. 

[13 type. Here we assume that the A (x) and B (x) are C and consider C 

J k J 

coordinate changes on IR n . 

Texts on partial differential equations develop the theory o f canonical 
forms for second order i i near partial differential equations i r* two var i 3 5 = 
(see Garabedian [33 and Courant anti Hilber t [33). Extensions to t l i: er -■ 
more variables are due to Cotton [^3 and Fredricks [53. These it.-uii.s 

CO 

involve finding C coordinates in which the principal part coefficients u ( 
become constant. This is always possible in two dimensions wi in suitable 


necessary and sufficient conditions in three or more dimensions, 


Little 


attention is paid to the first order coefficients B ( x ) after the new 

j 

coordinates are introduced. 

If there are C coor d 1 nates under which (i) becomes constant coefficient 

n .2 n 


< 2 > 


du 


+ v a _ ^ 

dt y . jk dx dx, , “i 

j.k-1 J j k j=l 


.2 

a u <*j , . 

+ z b - ~ f ( x , t > 

dx 


then Fourier transforms in the spatial variables yield 
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an ordinary differential equation in LI. 


Similarly, if there are C coordinates in which <1) takes the form 


- n .2 n 

au a u , oa 

dt . ^ a jk dx ,£x f . jk”j dx~ 
j>k=l j k j , k = 1 ' k 


f<x,t) , 


where a . and b , are constants, then Hormander [61 shows that spatial 
jk jk 


l-ourier transforms lead to 


Oil r*U 

— - 2 a , £ t, U - 2 b £ -r-— = Mt.t), 

c't . . jk-j-k . . , jk *k 


j * k = 1 


Under generic conditions, this first order linear partial differential 


equation can be solved bv the method of char 


( i . e . ordmf 


differential equations) and inverse Fourier transforms. The most important 


case involves a positive semidefimte matrix 'a. . ; of constant re: 


shall call the corr espond i nq equation a Kolmogorov equation, Hormance. 


writes the spatial partial differential operator in (h) as 


X. and X „ are C vector fields on IF . in fact, if the (A., 1 
i U ,1 k 


U) is symmetric, constant rank m, and positive semidefimte, the spat 


operator can be written as 


X + X _ . 
1 u 


For general partial differential operators of the form (6) Hormander proceed- 


to prove necessary and sufficient conditions for hvpoe 1 1 i p t i c i ty (i.e. C 


right hand sides imply C solutions). We require that our Kolmogorov 


equations be hypoe 1 1 i p t i c . Weber [73 constructed fundamental solutions for a 


class of equations related to those in (A). 


In this light the problems considered in this paper are: 


i) Given the partial differential operator i6) find necessary and 
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CO 

sufficient conditions so that there exist nonsingular C coordinate 
changes (local-near the origin) on IJ<P under which (6) becomes a 
constant coefficient partial differential operator. Standard 
differential geometry results (e.g. Spivak C83) are employed and 
the results are of no surprise. The conditions are derived here 
for the sake of completeness. 

ii) Given the partial differential operator (6) find necessary and 

CO 

sufficient conditions so that there exist nonsingular C coordinate 
changes (local-near the origin) on IR n under which X.X_,...,X in 
(6) are transformed to constant vector fields and X^ becomes a 
linear vector field. This makes the partial differential operator 
(6) of Ko 1 mogorov type , if the hypoel 1 ipt ic l ty conditions of 
Hormander are satisfied. 

The spatial operator in ( ^ ) is a Kg 1 mogorov operator if ii - 
hypoel 1 i pt i c . 

We remark that both problems i) and ii) can be generalized u< fne 

m p . ... 

partial differential operators T X L - 2 X “ + X., but we shall concentrate 

• j=i 1 

on the form (6). Moreover , we assume that X, , , . . . * X are i i 

1 P m 

i ndependent . 

Our principle tools are taken from the field of systems and control. 
However, the purpose of this paper is not to draw a parallel between 
controllability of systems of nonlinear ordinary differential equations and 
hypoellipl icity of partial differential equations, as this has been well 
established in the literature and in conference presentations. 

As we mentioned previously, problem i) is straightforward. Our work on 
problem ii) is analogous to the study of coordinate changes to transform a 
nonlinear control system on (R n , 
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(7) x = f<x)+Zug(x) 

j-1 j j 

to an n-dimensional controllable linear system 

< 8 ) x = F x + Gu . 

Here x = , f,g ,g , — ,g are C vector fields on R n , f<0) = 0, u = 

at id m 

) consists of real-valued functions? F is an nxn constant 

id m 

matrix, and G is an nxm constant matrix. Also f,F, and u are obviously 
different objects in our control discussion than in our p.d.e. discussion. 

We rely heavily on the results of Krener L9] and Responded [10], and, 1 n 
f ac t , our research essentially moves their results from the ordinary 

differential equation setting to the partial differential equation settir.c. 
Nonlinear control system (7) is replaced bv equation (6) with X tai-iuu ih« 
place of q , i = 1,2, . ..*m, and X.. taking the olar.n of f. The 1 i s vs ten 

i U 

(8) is replaced by the Kolmogorov partial differential operator 

m 2 c 

<9) t y f X , 

j=? 3 ° 

where each X is a constant vector -field and X,. is linear. Here X , X , . . ; , r. 

j 0 ; d :L 

correspond to the m columns of G and X^ corresponds Fx in (8). Wo want tht- 

span of the Lie brackets X , fX_,X 3 ad n * X . , X ), j=l,2,....m to be IF*’. 

j O j o i 

so we make the corresponding assumptions on the Lie brackets of vector fields 

in (6). As noted before we also suppose that X «X_,...»X are 1 inearl*. 

1 d m 

independent . 

Section 2 of this paper contains basic definitions and consideration of 
those linear partial differential operators, which can be made constant 
coefficient. In section 3, necessary and sufficient conditions are derived 
which classify those linear partial differential operators that can be moved 


to the Kolmogorov type. 
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II. Constant Coefficient Operators 

We begin with a set of appropriate definitions. 

If X and Y are C™ vector fields on IR n , then the Lie bracket of X and Y 


is 


CX,Y3 - | Y -X - Y * 
dx a>: 

where ~~ and are Jacobian matrices, x beinq the variable for IP . 

d x <?x 


Successive Lie brackets such as 


X, [X,Y] 


Y, C X , Y j 


X . [ X , Y 3 ] » Y 

J 


etc. can be taken. A standard notation is 


(ad X,Y) - \ 


<ad~ X,Y) = C X , Y J 


; ad a , Y ) ~ 


X , [ X , Y ] ; 


(ad J X,Y) = CX,(ad J 1 X , Y > 1 

We let <*,•> denote the dual product of one forms and vector fields. 
Given a C function h on R 1 we define the Lie derivative of h with respect tr 


the vector field X a? 


Successive Lie derivatives are 


L h = <rih , f > . 


L x h - h 


L h = L h 

A * 


L x h - L x L x h 


L x h ~ L x L x h ‘ 

Moreover, the Lie derivative of the one form dh with respect to X is 
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L x <dh) 


-e 


d(dh)* 


dx 


* 

+ (dh )-^— , 
dx 


where * denotes transpose. 

The three types of Lie derivatives satisfy the formula 


< 10) 


L <dh,Y> = <L (dh),Y> + <dh,CX,Y]>. 

A A 


We motivate our study by the following example. 


Examp 1 e g,l . Consider the partial differential operator 

.2 


1 1 ) 


* ^ ,, * 

(1 + x ; -- + (!+>• } + ?_>' + 

1 d ) ay. 3 dx c/x ... 

7X . 1 2 d 


1 


or IP: . The local coordinate change (near the origin) 


y = k n ( 1 + x i 
1 1 


( 12 ) 


Yq = 
c c 


v 3 = v : 


mo v e s (11) to t h e c o n s t- a n t c oe f f i c i e n t + o r m 


( 13) 


0 


, c b v 
<*y j 3 


This is discovered in the following way. First we write 


( 1 H ) 

and set 


n i I & n 

( 1 + x . ) - — (1+x, ) + 2x„ — — *• + 

1 a * j l 1 j I 3 r ''x dx 


( 15) 


X 1 = 


Vx i 


0 

0 

, X. = 

2x 


n 

3 

0 


1 


Thus (11) becomes 
(16) 


X, + X_ 
1 U 


Since and X^ are linearly independent and the Lie bracket C X ^ 


- 0* standard differential geometry results (see C 8 ] ) imply that 
transformat ion (12) takes 


X to 

1 

0 

and to 

0 

0 

1 

0 

w 

0 

1 


- 6 - 






] 

the 




Hence* the partial differential operator (11) is moved to the constant 
coefficient form (13). 


We now prove our result concerning transformation to constant 
coefficient operators. Again, this result is trivial from the differential 
geometry viewpoint . 

co n 
Theorem 2 . 1 . Given the C partial differentia] operator on IR 


( 17) 


Z X + X , 
j=i J 0 


where X. * X_ , „ . . X , are linearly independent for m<n and X . , . . . , X are 

0 1 2 m 1 c m 

linearly independent for m=n, there exists a non-singular (local) coordinate 

change on fR P under which X .X,j...,X become constant vector fields if and 

0 1 m 

o n 1 y if 


i i 8 J 


,x s! 


- 0 for all U<r,s-.m. 


Proof . If m<n results from [0] indicate there are nonsinqular coord i na m 
t a k i nq 


0 

0 


I i 


I 

r 

j y. 

i 

1 

s 




I 

} 

r, 


• 

I 

; . >• tO 

- 

w . - X to 

I * 

1 <■- 


m 

I 1 

| 

1 

j 


j 

0 
i o 

n 

I 

0 

.. " j 

X . to 
0 

1 

0 

, , . th 

— ( n-m- 1 ) pi ace 


th , 

< n-m 1 piece 


if and only if (10) is satisfied, 
If m = n, we can move 
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' o ’ 


' 0 


■ i 


0 


0 


0 

to 

* 

X to 

• 

i . . . > X to 

* 


- 

2 

. 

m 

• 


0 


1 


0 


1 


0 


0 


if and only if 


= 0 for all l<r»s<m. Setting 


we find the only possible nonzero column of 


<J7 


dx 


07, 


f. i i dx 

i ? . > X I i s j J i 

I 3 °! | ■ ! 


! . ! 

I 


I - - n i 
! ox i 
i. J J 


■for i “ 1 . . . »n. Then X. is a constant vector field if and onl v il X % >. • 

O j i '••• ; 

£ 0 , } = ] , d , . . . , n . ;j 

We now study our problem ii), the main consideration of this poDt?- . 


Ill- Kolmogorov Operators 

^ 2 

We examine the partial differential operator (6) v. X + X where 

,=i J 

X »X are linearly independent and X. vanishes at the oriqin in it' 1 ', 

i d m s' 

We derive conditions under which C coordinate changes (local -near the 
origin) exist taking (6) to the Kolmooorov operator (9). 

As stressed in the introduction, the main contribution of this paper is 
realizing that the results of Krener L9J and Respondek CIO] in the nonlineai 
systems (o.d.e.) and control area can be applied to partial differential 
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operators 

< 19) 


n £ n a 

Z A . (X) + ZB < X ) rjL- 

jk dx .dx. . , j dx . 

J>*<=1 J k J = 1 J 


The linear controllable system 

(20) x = Fx + Gu 

is replaced by the hypoelliptic Kolmogorov partial differential operator 

(with a , and b , constant) 

.1 k jk 


(21 ) 


n .2 n 

n 

2 a , — — + Z b , >: ~z . 

, , j k ax ax, . , , jk j ax, 

j*k = J j k jtk=l ‘ J k 


Of course we shall study operators (19) and (21) in our vector field 
notation. First we wish to examine parallels between system ‘.20; and 
operator (21 ) . 

For the control system (20) the Kronecker indices and eigenvalues ot the 
F matrices are invariants under coordinate changes. For the operator (21/ re- 
introduce Kolmogorov indices and note that these and the eigenvalues of the L* 
matrix are invariants. Canonical forms which parallel controllable canonical 
forms for (20) will occur in our work. 

If the matrix A - (a , ) in (21) has rank m, let X. , X . X_ ♦ . . . > X be 

j k 0 1 d m 

m _2 

linearly independent vector fields so that (21) becomes Z X. + X_. We set B 


3=1 


j 


= Jacobian matrix of the vector field X_. If the partial differential 

0 

operator in (21) is hypoelliptic (in this case X >BX S' , n = 

} j j 

l,2,...»m, span IP ) we introduce the following process: 
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1) Write out the grid 

X. 


”1 

A £ • • • 

m 

BX 'l 

bx 2 . . . 

BX 

m 


6% • • • 

B 2 X 

m 


B"' 1 X. 


b" . . . b" ‘x 

£ m 


Start at X^ and move left to right across the first row, then start 


at BX and move across the second row, etc 


Throw out any vector field that is a linear combination ot the 
preceding vector fields in the grid. Discard all vector fields r» . 
the column below this vector - field. 


Continue until n linearly independent vector fields die found and 
all others have been discarded. 


Let <- = number of entries remaining in the ^ " column of the arid- 

J 


j - 1 >£*..., m . 


6) Renumber X,,X % if necessary? so that V f ' . 

1 £ fit i “ q " “ m 


Def . The integers 9 are called the Kolmoqorov indices of 

- 1£ m a 


partial differential operator (£1). If the vector field notation for i£l) i- 
m £ 


<9; Z X + X. we also have the inteqers (\ f associated with (9). 

, i 0 * 3 £ m 

j = 1 ‘ 


By assumption + €_ +...+( = n. 

Id m 


We ask if the partial differential operator (6) Z X7 4 X is a 


coordinate change away from the Kolmogorov partial differential operator (9) 
m 


Z X. + X_ having given indices . Before stating the qeneral 

js=i j 0 1 £ m 




theorem we present an example. 

Example 3.1 . Consider the partial differential operator 


d 2 


(22) *tx — — + *tx ^ ^ — , — — 

3 _ 2 3 ax_ax n ,2 

dx n 2 3 <?x_ 

d 3 


^ + ( 2+x_ ) ■ ^ ■ + fx_ - x^ + 2 ( x _-x^~ ) x 


3 ax. 


2 3 3 ax. 


on 


We write this as 


(23) 

Let t i nq 


n a d 
2 x_ + 


„ a a 

2x_ - + - — I + 

3 x dx J 


0 


a 

' 

OJ 

2x 3 

5 X — 

x 2 -x 3+ 2(x 2 -x 3 )x 3 

X _ 

1 

0 

3 



0 


3 ax ax 

2 3- 


x i = 


we find (23) becomes 

<2h; Xj - X 0 . 

The local coordinate chance' 


3 _ 2 h 

*1 * v*2 * e -a y 3 ' y 3 


X _ = X - X _ = L X 

2 2 3 X„ ! 

0 


* 3 = X 3 = L X . X 2 


*2 " *3 + 2U S -4 ] *3 


— — + x„ — — 
c»x 3 dx 


take (24) to 
(25) 


V * *0 ' 
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where 

X, = 

~ 0 ! 
0 

j 

and X ^ = 

x_ 

d 

>: _ 


1 

1 

0 

3 

0 j 


This yields the Kolmogorov partial differential operator 


(26) 


d -d-d 

- 8 8 .- 3 

dx 3 dxj dx g 


Since 
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B = 


0 1 0 
0 0 1 
0 0 0 


2- 3 

Xj, BX^, B X^ span [R , and the single Kolmogorov index is ^ = 3 = n. 
In equation (B 4 *) we remark that 


and 


X .CX ,X 3, <ad 2 X .X > span (R 3 (near ( 0 , 0 , 0 )) 
C < ad r X , X ) , ( ad S X Q , X j ) 3 s for CKr,s<3 


We now state and prove our main result. The general proof will follow 
from analogous results from C9] and [ 1 0 J » but we shall present a proof in 
the case m = 1 for some sort of completeness. 

m 

cl 

Theorem 3.1 . The linear partial differential operator (6) Z X + with 

j-J 

m 

X, 9 X_ > . . • * X linearly independent on fp‘ ' and X.(0) = 0, can be transformed bv 
1 dm 0 

nonsingular coordinate changes (local-near the origin) to the Kolmogorov 

m - 2 

partial differential operator (9) Z X + X_ havmq indices * » . . . » t if 

, i 0 1 B m 

j-1 


and only if 


J r i 1 

the set { X ,[ X_ , (ad X ,X), X_, [ X. « X„:i , , . . 

\ 1 0 1 0 1 B 0 B 

fp-1 € -1 

( ad X * X )>...> X » [ X pi ( ad X > X ) > 

OB m 0 m 0 ml 


is 1 i near 1 y i ndependent . 

and m) the Lie brackets of every pan of vectors fields in 
j' , *1 # 

\ X » [ X * a ] * . . . > ( a d X ? X ) > X 9 [ X ♦ j * . . . » ( ad X » X _ ) < 
1 O 1 0 1 B 0 d O d 


X >[X ,X (ad X i X > ) is zero, 

m 0 m 0 m 


Proof. (For m = 1» and the operator X, + ) ; 

1 3 0 

Since X ^ » [ X ^ 9 X ] 9 - - - > ( ad n *X^*X ) are linearly independent and have zero 


* 12 - 


Lie brackets* we have coordinates so that 


original page is 

0F p °OR quality 


Xj becomes 


0 

0 

0 

1 


[ X * X ] becomes [X.*X.3 
0 1 0 1 


0 

0 

1 

0 


If x o = 


(ad n becomes (ad n *X 


n -1 

'n 


i 

0 


0 

0 


si <\ 

* then CX ,X 3 = 0 implies — -- = = 

fi x 

n 


i n~c * n * 


(28) 


2 - - ^ 

« ad X rj 9 X ) = 0 implies = u , i - 1 > 2 , . . . % n -3 , n~ 1 > n , - 

C*k dr 

n- i 


n- 1 - - jK 'i U,J ] 

(ad X ,X > = 0 implies ----- — - Os i = ^ , 3 . . . n , 

<. r >, ; .f >' 

c r' 


<7 o . 

Now C(ad n X *X )»X 1 = 0 yields — 

Oil - 

o' X .-7 >. 

i n 


0 % l » i - 1 i 2 ... s n 


d 2 & 


L(ad X Q ,X 1 ) 9 [X^»Xj31 - 0 yields 


- — 0, i » j - 1 » 2 * . . . » n 


dx dx 
j n -1 


( 29 ) 
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C(ad n X 0 ,X 1 ), (ad" *X ,X )3 = 0 yields — - — -i 

dx dx. 
j i 


~ 0 * i » j — l,2,...,n . 


- 2 


Hence is a linear vector field as promised and is a Kolmogorov 

partial differential operator. a 

In the above proof 


(30) 


0 

0 

0 

1 


and X, 


- 2 


11 

D 21 


(n-1 ) 1 
3 ni 


1 0 ... 0 

0 1 ... 0 

> 0 ... 1 

0 . . . 0 


and + X^ is in a canonical form. If we set 



f X 1 

i 


x ? 

! 

. 


Vi 


x I 

n j 


: i = x i 


2 ' L X . ' i 
u 


r “ L. - , 

n X .. n- i 

i v > 

we have a canonical form m which the matrix defined h v the new X. is i. < 

(.) 

rational canonical form and is as before. If nr: . 1 , we get the analogue of 
c g n t r o 1 1 ab 1 e canonical form. 

Ue have developed a theory giving necessary and sufficient conditions 
that a second order linear partial differential operator he a coordinate 
change away from a Kolmogorov operator. 

Future research will be in two directions: 

1) Expand the transformat ions used to include ’’appropriate types" of 
feedback. This research is presently underway, and first thoughts 
were to include results in this paper. However , the process of 


feedback , 


as applied in transformation theory, is not well 


addressed in the partial differential equation literature* where 
coordinate changes are standard fare. Therefore* we decided that 
separate papers are appropriate. 

2) Extend all results to the discrete setting. A Ph.D. student of the 
first author is currently working on this project. 
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